Water vapor at a translational temperature of one kelvin. 
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We report the creation of a confined slow beam of heavy- water (D2O) molecules with a transla- 
tional temperature around I kelvin. This is achieved by filtering slow D2O from a thermal ensemble 
with inhomogeneous static electric fields exploiting the quadratic Stark shift of D 2 0. All previous 
demonstrations of electric field manipulation of cold dipolar molecules rely on a predominantly lin- 
ear Stark shift. Further, on the basis of elementary molecular properties and our filtering technique 
we argue that our D2O beam contains molecules in only a few ro- vibrational states. 
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Cold dilute molecular systems are rapidly emerging as 
a front line area at the interface of quantum optics and 
condensed matter physics 0]. An increasing subset of 
this activity centers around the creation of cold dilute 
gases of molecules possessing electric dipole moments. 
These in particular, owing to their long-range anisotropic 
interaction, hold the promise of novel physics, where two- 
and many-body quantum properties can be systemati- 
cally studied. Cold dilute gases of dipolar molecules can 
be produced by forging a tight bond between two chemi- 
cally distinct species of laser-cooled atoms, e.g. RbCs 0- 
Alternatively, cold dilute gas ensembles can be created 
by buffer-gas loading Q or electric-field manipulation of 
naturally occurring molecules like ND 3 0,0, H 2 CO 3, 
metastables like CO Q or radicals like YbF ||, OH [<j, 
NH 01 ■ So far all the cold molecules made available 
with electric-field-based methods have a Stark effect (in 
their relevant states) which is predominantly linear in the 
important range up to 150kV/cm. 

Here we report the creation of a slow beam of heavy- 
water (D2O) molecules, which experience a quadratic 
Stark effect. The cold D2O molecules are filtered from 
a room-temperature thermal gas and have a transla- 
tional temperature around 1 kelvin. Because the Stark 
shifts are quadratic in the electric field, it follows that 
forces exerted by inhomogeneous electric fields are rel- 
atively small for D2O compared to molecules with sim- 
ilar dipole moments but with linear Stark shifts. R is 
therefore by no means obvious that significant quanti- 
ties of slow D 2 molecules can be produced by means 
of electric-field-based methods. Our experimental re- 
sult therefore underlines the versatility of the velocity- 
filtering method. R is an enabling step towards future 
trapping of molecules for which the ratio of elastic to in- 
elastic collisions is expected to be more favorable than for 
molecules with linear Stark shifts 0] . An additional ad- 
vantage of the quadratically Stark-shifted molecules like 



D2O is the possibility to perform precise spectroscopic 
measurements insensitive to stray electric fields, to the 
first order. Moreover, water is abundant in interstellar 
space at low densities and temperatures from a few kelvin 
upward, playing an important role in the chemistry of 
molecular clouds 0| . The conditions in these clouds are 
remarkably close to those achieved in our experiment, 
opening up the possibility to investigate in the labora- 
tory chemical reactions under conditions found in space. 

This Letter is structured as follows. First we discuss 
general features of Stark shifts of molecular states with 
particular references to D2O. We then present our exper- 
imental work with D2O. This is followed by arguing from 
first principles that the resulting beam of D2O is domi- 
nated by only 4 rotational states, despite starting with a 
thermal source of molecules at 300 K. 

Several techniques have recently been developed to ma- 
nipulate molecules Q with electric fields. All of these ex- 
ploit the Stark effect to exert a force on the molecules. In 
contrast to atoms, molecules can have a permanent elec- 
tric dipole moment. Such molecules have much larger 
Stark shifts than non-polar molecules. However, as de- 
scribed below, a large dipole moment alone is not enough 
to have a strong Stark effect. The direction and magni- 
tude of the force exerted on the molecule in an inhomo- 
geneous electric field depends on the details of the molec- 
ular rotational state. Assuming the Stark shift to be a 
monotonic function of the electric field, the molecule can 
be either in low-field-seeking (LFS) or high-field-seeking 
(HFS) states, depending on the sign of the Stark shift. 

The condition for having a linear Stark effect is that 
the component of the dipole moment d along a space- 
fixed direction, we can choose z, is non-vanishing (i.e. 
(d ■ z) ^ 0). Strictly speaking this requires a finite elec- 
tric field, but it can be arbitrarily low. The linear Stark 
shift is typically found in symmetric top molecules and 
is proportional to \d\KM, with K and M representing 
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the projection of the total angular momentum J on the 
molecular symmetry axis and on the z-axis, respectively. 
Of course, no first-order Stark effect occurs if either M 
or K - or both - are zero. If the degeneracy in zero 
field is lifted, e.g. by fine-structure splitting, inversion 
doubling, nuclear quadrupole interaction in a symmetric- 
top molecule or A-doubling in linear molecules, the Stark 
splitting will no longer be linear in the limit of zero field. 
However, often those interactions are small enough as 
to lead to a nearly linear Stark splitting in the applied 
electric field range. 

In general, the linear Stark effect condition, (d- z) ^ 0, 
is not fulfilled in asymmetric top molecules. Under cer- 
tain conditions, however, polar asymmetric molecules can 
also exhibit (nearly) linear Stark shifts. If the asymme- 
try is weak, the states that correspond to K ^ in the 
prolate or oblate top limit, will always be close to being 
degenerate. Those states show nearly linear Stark shifts 
if they are coupled by the Stark interaction. This is the 
case if the dipole is along the a-axis in the prolate limit, 
or along the c-axis in the oblate limit, where we follow the 
convention to label the axis with the smallest moment of 
inertia and hence the largest rotational constant with a, 
the intermediate axis b and the axis with the largest mo- 
ment of inertia c. An example is the nearly symmetric 
(prolate) top molecule H2CO. 
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FIG. 1: The lowest rotational energies of D2O as a function 
of the applied electric field E. The 5 most abundant states 
I J, r, M > in the guided beam are indicated. The Stark shifts 
are obtained by numerically diagonalizing the Stark Hamilto- 
nian (J = ... 12), following Ref. Q. 



True asymmetric top molecules in general have 
quadratic Stark shifts. Exceptions can occur for some 
states, if the dipole is oriented along the axis of largest 
or smallest moment of inertia. For molecules with their 
dipole oriented along the &-axis, we found no exceptions: 
all rotational states have a non- linear Stark shift. Wa- 
ter, both H 2 and D 2 0, presents such a case and a few 
levels of D 2 are depicted in Fig.Q The quadratic be- 
havior is obvious; only for the highest most abundant 
states, the | J = 3, r, M) states (where r is a pseudo quan- 



tum number labelling the state), a deviation is found. 
Moreover, the large rotational constants ^ij f° r D2O, 
A = 15.394 cm" 1 , B = 7.2630 cm" 1 , C = 4.8520 cm" 1 , 
imply large rotational level spacings. Avoided level cross- 
ings are neither expected nor found, and second-order 
perturbation theory is a reasonable approximation for 
the Stark shift computation of H2O and D2O. Moreover, 
since the contribution to the perturbation from each cou- 
pled pair of states is inversely proportional to the energy 
gap between the pair, the shift will be proportional to the 
density of (rotational) states, and therefore very small 
for the sparse rotational spectrum of D2O and H2O. Our 
choice of working with D2O as opposed to H2O has partly 
to do with the larger Stark shifts of D2O because of its 
smaller rotational constants. The treatment of the gen- 
eral case of an asymmetric molecule, where the dipole 
is not necessarily along one of the principal axes, is, of 
course, more involved. 




FIG. 2: (Color online) Schematic of the experiment. On 
the left is the effusive source, which injects thermal D2O 
molecules into the 4-wire guide. Neighboring electrodes 
have opposite polarity, creating a quadrupolar electric field. 
Molecules that are slow enough are guided through the first 
and second (not shown) 90° bends and are finally detected by 
a mass spectrometer. 

Our apparatus is depicted in Fig. It consists of 
a room-temperature effusive thermal source, which in- 
jects D2O molecules directly between four 50 cm long 
electrodes set up in a quadrupole arrangement, with 
neighboring electrodes having opposite polarities. The 
guide has two 90° bends with a radius of curvature of 
25 mm. The quadrupolar electric field defines a two- 
dimensional potential well. This well has a depth that 
depends on the internal molecular state, e.g. for the 
\J,t,M) — |1,1,1) state with a positive Stark shift of 
0.20 cm" 1 at lOOkV/cm, the depth amounts to 0.29 K. 
Molecules with transverse kinetic energy exceeding the 
potential depth escape the guide. In the bends the lon- 
gitudinally fast molecules escape while the slow ones are 
kept due to the action of the centripetal force. These 
are guided through two differential pumping regions into 
an ultrahigh vacuum chamber for mass-spectrometric de- 
tection at the end of the electrodes. Heavy water is con- 
venient for this purpose, as the background at its mass 
is virtually zero. The longitudinal velocity distribution 
of the guided D2O beam was determined by a time-of- 
flight method [|| at an escape field, E, of 115kV/cm. We 
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found a most-probable velocity of 24m/s in the labora- 
tory frame, corresponding to a longitudinal temperature 
of as 1.4 K. The transverse temperature is expected to 
be on the order of 0.1 K, as the guide presents a smaller 
transverse velocity cutoff value than the corresponding 
longitudinal velocity cutoff. 

The flux dependence on the escape field, E, and 
hence on the applied electrode voltage V, is character- 
istic of the nature of the guided molecules' Stark shift. 
This can be seen as follows: Let x,y be directions or- 
thogonal to and z be parallel to the quadrupolar axis. 
Let f Vxy , be functions proportional to the flux cross- 
ing the planes of unit area perpendicular to the x, y, z 
axes, respectively 0. Then f Vxy oc exp(— v 2 y /a 2 ) 
is bi-directional and f Vz oc v z cxp(— v 2 /a 2 ) is uni- 
directional along the positive z axis. Here, a = 
^j2ksT /m, /cb the Boltzmann constant, T the temper- 
ature of the reservoir where the beam originates from 
and m the molecular mass. Hence the total guided flux 
$ oc Jo*'™ 1 dv x j^ v ' m ^ dv y Jo i,max dv z f Vm f Vy f Vz , where 
Vx,y,z,max are the maximal guided velocities in each direc- 
tion. As a > V XtVlZtmaK , $ OC Vx,max%max«z, max - TllC 

maximum kinetic energy, t/k,max is given by the escape 
energy of the guide, i.e., the Stark shift in E, which is 
proportional to the applied electrode voltage, V. Hence 
for molecules with a linear Stark shift, <& oc U 2 oc V 2 . 
For molecules with a quadratic Stark shift, f/k oc V 2 and 
hence, $ oc U 2 oc V 4 . 

Our detector, a quadrupole mass spectrometer, con- 
verts molecules to ions by electron-impact ionization. 
The ions are then mass-selected. Measurement on vari- 
ous molecules with a linear Stark shift indicate that the 
signal is to a good approximation proportional to the 
guided flux 0, Q . Scaling (for detector counting efficien- 
cies, angular divergence of the beam exiting the guide) 
and corrections (velocity dependent detection, branch- 
ing ratios of ionization) are needed to convert our mea- 
sured count rates (plotted in Fig. to the absolute flux 
as 7 x 10 7 s _1 at V = 7kV, corresponding to an electric 
field depth of the guide of E = 134kV/cm. The error 
margin of the flux is estimated to be of the order of a 
factor 2. The quartic dependence on E is clearly visi- 
ble in Fig. |21 and proofs the quadratic Stark shift of the 
guided molecules. Indeed, with the same apparatus it 
has been observed 0,0 that for H 2 CO and ND 3 (linear 
Stark molecules), the flux depends quadratically on V. 

As our slow beam originates from a room-temperature 
source, many rotational states are populated. This is 
illustrated in Fig. 01 where the Stark shifts of D2O in 
a field of lOOkV/cm have been plotted as a function of 
the zero-field rotational energy. The Stark shifts are ob- 
tained by numerically diagonalizing the Stark Hamilto- 
nian for a rigid asymmetric rotor, following the procedure 
of Ref. . It is known that the rigid-rotor assumption 
is only a coarse approximation when estimating the abso- 
lute energies of D2O states. In the present case, however, 
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FIG. 3: Detector signal versus the electrode voltage V. The 
data follow a quartic law in V, as illustrated by the V 4 fit 
(solid line). The dotted line shows a V 2 fit attempt. 



the approximation is expected to be good if one is only 
interested in the Stark shifts and not in the absolute ener- 
gies, because the Stark shifts are caused by the coupling 
of adjacent states that do not differ much in their sets of 
rotational quantum numbers, leading to relatively small 
sensitivity to the centrifugal distortion. We have also ne- 
glected hyperfine couplings, which is completely justified 
in the range of field strengths used in our experiments. 
As input for the Stark shift calculations we took the rota- 
tional constants and the dipole moment fj, = I.87Debye, 
which is directed along the 6-axis Note that as a 

general trend, the Stark shifts decrease with rotational 
energy. 
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FIG. 4: Stark shifts, AWs, in an electric field of lOOkV/cm of 
LFS rotational states of D2O in its vibrational and electronic 
ground state, as a function of the zero-field rotational energy, 
E TO t- The inset shows the Boltzmann factor at T=300K. The 
five most Stark-shifted states are labeled | J, r, M) . Note that 
the spin statistical weighting is not shown. 

As discussed earlier the room temperature rotational 
spectrum of D2O is very sparse. The populations in the 
thermal gas in the effusive source can be estimated by 
weighting all thermally populated states with the Boltz- 
mann factor, their orientational (M) degeneracy and the 
spin-statistical weight |l6j . Further, only very few of 
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State Contribution [ % ] E Iot [cm -1 ] AWs [cm -1 ] 
|3, —2, 1) 21 74.53 0.16 

|1, 1, 1) 21 22.66 0.20 

|2, 0,2) 17 49.30 0.13 

|3,-2,0) 13 74.53 0.18 

|3,-2,2) 8 74.53 0.10 

TABLE I: The most dominant rotational states \J,t,M) of 
D2O with their partial contribution to the flux, their zero-field 
energy and the Stark shift AWs at E = lOOkV/cm. 



these levels have large enough Stark shifts to be guided. 
This selection is much more pronounced for molecules 
with a quadratic Stark shift than for molecules with a lin- 
ear Stark shift. The intuitive reason is that for molecules 
with a quadratic Stark shift, the electric field must first 
orient the dipole in space, which is harder for faster ro- 
tatin g m olecules. Indeed, for these molecules the Stark 
shift [lj| is approximately proportional to 1/(J + 1). 
Hence, the (maximum) Stark shifts decreases with J, as 
can be seen for D 2 in Fig. 0] This should be compared 
with the dependence of the Stark shift on J of molecules 
with a linear Stark shift. For the generic example of a 
symmetric top, this shift is oc \d\KME/(J(J + 1)) 0. 
As (K,M) e {—J, . . . , 0, . . . , J}, the Stark shift of the 
maximum (K, M) for molecules with a linear Stark shift 
will not decrease with J. In fact, knowing that the flux of 
D 2 molecules is proportional to the square of the Stark 
shift, and assuming that this dependence holds for each 
state, the four most populated states contribute more 
than 70% of our guided flux. The partial contributions, 
zero-field energy and Stark shifts of the 5 most abundant 
states are summarized in table [I] One should note that 
the beam purity [TtT ] is independent of voltage changes 
as long as all the states are quadratic in nature. 

In conclusion, we have demonstrated the effective 
Stark manipulation of a polar molecule with quadratic 
Stark shifts over the range of applied fields of 0- 
135kV/cm. This experimentally shows the feasibility 
of the velocity-filtering method for quadratically Stark- 
shifted molecular states. Using this method we have cre- 
ated water vapor (D2O) at a translational temperature 
of « 1 kclvin. Its quadratic Stark effect combined with a 
large rotational spacing make D 2 a promising molecule 
for electric trapping [13 and even evaporative cooling. 
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